New US Patent Application 



Applicant: Netzsch-Geraetebau GmbH 

Title: Device for detecting thermal conductivity by 

of optical pulses 

Inventor: Dr. Opfermann and Juergen Blumm 



1 



Description 

Field of the Invention 

The present invention relates to a device for detecting thermal con- 
ductivity by optical pulse techniques. 

5 Prior Art 

Various devices for determining the thermal conductivity by optical 
pulse techniques have become known. For example, the German Pat- 
ent DE 33 17 513 C2 describes a device of the claimed general type, 
wherein energy is supplied to a sample by laser pulses. The infrared 

10 radiation emitted by the sample is measured by means of a pyrometer. 
To maintain those effects at a constant level, which are produced by the 
laser pulses on the measuring result, a time-limiting camera release is 
disposed between the laser and the sample. This camera release en- 
sures at least a reproducible pulse shape. This device entails the disad- 

15 vantage, however, that in the case of long pulses or extremely thin and 
thermally conductive samples, the results of measurement are falsified 
by the pulse shape of the laser. 

The US Patent 4,928,254 discloses a device for and a method of de- 
20 tecting the thermal conductivity by application of an optical pulse 
method. There, a sample is irradiated with a succession of optical 
pulses from a YAG laser. A sensor is used to detect the infrared radia- 
tion emitted by the sample. Then, the thermal conductivity is calculated 
by deriving it from the result of measurement by means of a computer. 
25 Moreover, a certain fraction is decoupled from the laser beam and sup- 
plied to a photo diode. The computer can then determine the number of 
laser pulses issued on the basis of the signal from the photo diode. In 
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this method, too, the shape of the laser pulses it left out of considera- 
tion. 

The European standard EN 821-2 establishes a maximum period of a 
5 laser pulse at 1 % of the half-rise time (the period elapsing until 50 % of 
the maximum signal level is reached) of the detection signal of the infra- 
red sensor in order to avoid a high impairment of the result of meas- 
urement beyond proportion, which is caused by the influences taken by 
the pulse lengths of the laser pulses. 

10 

In order to permit also measurements with longer laser pulses or in 
cases of a higher thermal conductivity, the pertinent literature proposes 
the approximation of the laser pulse and the correction of the detected 
curve shape of the infrared sensor by performing a convolution with the 

15 approximated laser pulse. A particularly expedient approximation tech- 
nique is proposed by Larson and Koyama in: Journal of Applied Phys- 
ics, 38 (1967), p. 465. According to that publication, an exponential form 
of the curve shape of the laser pulse is assumed. With this approxima- 
tion in correspondence with the European Standard EN 821-2, it is pos- 

20 sible to realize a maximum period of a laser pulse, which corresponds 
to 10% of the half-rise time of the measuring signal of the infrared sen- 
sor. Measurement in a range going therebeyond is not possible at pre- 
sent. 

25 

Brief Description of the Invention 

The present invention is based on the problem of improving a device for 
or a method of detecting the thermal conductivity by means of optical 
pulse techniques in such a way that even in the case of comparatively 
3 o long optical pulses as well as extremely thin samples or of an extremely 
high thermal conductivity, measurements can be performed. This re- 
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lates particularly the range in which the maximum period of a laser 
pulse exceeds 10 % of the rise time of the measuring signal of the infra- 
red sensor. 

5 In accordance with the present invention, this problem is solved by the 
provisions defined in the independent Claims. Expedient improvements 
of the invention are the subject matters of the dependent further Claims. 

An inventive device for detecting the thermal conductivity of a sample 

10 (1 ) comprises a furnace (2) for heating the sample to a previously se- 
lected temperature. Moreover, a source of radiation or light source (5) is 
provided, preferably in the form of a laser or a flashlight for emitting ra- 
diation (3) for heating the sample. It is equally possible to emit also any 
other kind of electromagnetic radiation instead of light. The radiation is 

is preferably emitted in the form of short pulses. The energy absorbed in 
the sample results in heating of the sample. In correspondence 
therewith, the infrared radiation (4), which is emitted by the sample, var- 
ies, too. The time history of this infrared radiation is detected by means 
of an infrared sensor (6). The signal from the infrared sensor (6) is sup- 

20 plied, via a signal amplifier (12), to an analyzer unit (13). This analyzer 
unit is so designed that is determines the thermal conductivity of the 
sample by deriving it from the time history of the sensor signal. For ad- 
justment of the result of measurement by elimination of falsifications 
caused by the time history of the optical pulse, a device is provided for 

25 measuring the time history of the optical pulse. This device comprises a 
decoupling element (7) such as a semi-transparent mirror that decoup- 
les a predetermined fraction of the radiation of the light source (5) as a 
reference radiation and supplies the latter to a sensor (9) for measuring 
the laser radiation. This value so measured is then supplied to the ana- 

3 o lyzer unit (1 3) for detecting the time history of the optical pulse. In order 
to achieve a sufficient resolution of the sensor it is necessary that the 
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sensor has a bandwidth in correspondence with the invention, which is 
substantially wider than the reciprocal value of the pulse length of the 
radiation pulse. The bandwidth wider by more than a factor of 10 than 
the reciprocal value of the pulse length is preferred. The time history of 
5 the radiation pulse is expediently detected by digitizing and storing the 
measured values. The duty cycle must be selected here in such a man- 
ner that a sufficient number of measured values is determined, prefera- 
bly at least ten values, throughout the period of the radiation pulse. The 
analyzer unit is so designed for determining the thermal conductivity 
10 that it corrects initially the signal measured by the infrared sensor (6) by 
application of the signal of the sensor (9) for measuring the laser radia- 
tion. 

Tests have shown that with such a correction the ratio of the pulse 
15 length of the laser pulse to the half-rise time of the measuring signal 
can be extended up to a value of 0.75. Consequently, the measuring 
range can be widened by a factor higher than 7, compared against a 
correction in compliance with the standard EN 821-2. 

20 With such an inventive correction is has now become possible for the 
first time to measure multi-layer samples with a high precision. 

It is, of course, also possible to provide a device for cooling the sample, 
instead of the furnace (1 ). For simplification of the representation, the 
25 term "furnace" is referred to in this document because in the majority of 
applications such a furnace is actually provided. Such a furnace 
serves, however, the function of maintaining the temperature in the en- 
vironment of the sample at a constant predetermined value. 

3 o Furthermore, the analyzer unit (1 3) is designed for correcting the signal 
measured by the infrared sensor (6) by performing a convolution with 
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the signal of the sensor (9) for measuring the laser radiation. This pro- 
vision preferably serves to achieve the normalization to a corresponding 
comparative signal that relates to an excitation by an ideal laser pulse in 
the form of a Dirac pulse. 

5 

In addition, the analyzer unit (13) presents a configuration for approxi- 
mation of the laser pulse by sections. In this technique, the laser pulse 
is approximated in at least two sections by means of exponential func- 
tions. With typical laser pulses, the time history can be subdivided into 
10 the following sections: 

(a) a high-speed rise starting section 

(b) a generally slightly downward sloping plateau section, and 

(c) a high-speed decay section following the cut-off of the pumping light. 
In correspondence with the curve shape or the duration of the laser 

15 pulse relative to the rise times and the decay times, this pulse can be 
advantageously approximated by defining the sections (a), (b), (c) or, in 
the event of short laser pulses, only defining the sections (a) and (b). As 
a matter of course, in the case of particularly complex pulse shapes, it 
is also possible to approximate a laser pulse also by the subdivision into 

20 further sections. 

According to the invention, the analyzer unit (13) is designed for deter- 
mining the difference in Delay between the rated zero point in time and 
the starting time of the optical pulse. 

25 

Another embodiment of the analyzer unit (13) consists in the provision 
that for approximation of the intensity of the laser pulse /,(*) during the 

rapidly rising edge the following formula is used as a function of the 
time t: 

3 0 I x (t) = A • [l - exp {-(/ - Delay) I r, }] 
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wherein: 

A scaling factor 

Delay : time delay of the starting point of the upward edge rela- 
tive to a zero point in time 

r, : time constant of the rising ramp 
The variables of this function and the functions described below are de- 
termined by derivation from the measured signals of the signal of the 
sensor (9) for measuring the laser radiation. In general, good results are 
achieved already with this single-term formula. 

As a matter of fact, an improvement of the results can be achieved by 
application of the following formula for I { (t) : 



/,(/) = A 



1 — exp {-(* - Delay) I r I2 } + r " ■ exp {-(; - Delay) lr n ) 



r l7 -r, 

wherein 

is Delay : time delay of the starting point of the upward edge rela- 
tive to a zero point in time 
t u : first time constant of the rising ramp 
t u second time constant of the rising ramp 

20 A further configuration of the analyzer unit (13) consists in the provision 
that for approximation of the intensity of the laser pulse I 2 (t) during the 
generally slightly downward sloping plateau the following formula is ap- 
plied as a function of time t: 

I 2 (t) = /, (0 • exp {-(/ - Delay) I r 2 } 
25 wherein 

/, (0 : formula for defining the intensity of the laser pulse during 
the rising edge 
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Delay : time delay of the starting point of the rising edge relative 

to a zero point in time 
r 2 time constant of the rising ramp 

7,(0 is an optional term defining the previous time history of the inten- 
5 sity of the laser pulse. Preferably, one of the above-presented formulae 
is applied here. 

Another embodiment of the analyzer unit (13) consists in the provision 
that for approximation of the intensity of the laser pulse 7 3 (/) during the 
10 high-speed downward edge, after cut-off of the pumping light, the fol- 
lowing formula is applied as a function of time t: 

h (0 = i 2 if = Dela y + ) • ex P Hf - Dela y 

wherein 

I 2 (t) : formula defining the previous time history of the intensity 
15 of the laser pulse 

Delay : time delay of the starting point of the upward edge rela- 
tive to a zero point in time 
t e : difference in time between the beginning of the down- 
ward edge and the beginning of the upward edge 
20 r 3 : time constant of the downward edge. 

I 2 (t) is an optional term defining the previous time history of the inten- 
sity of the laser pulse. Here, one of the formulae presented above 
should be applied. 
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An improvement of the results can be achieved by the application of the 
following formula for I 3 (t) : 



h(t) = I 2 (t = Delay + t e ) 

• — 2s exp{-(/ - De/ay - /.)/r 32 } ^ exp{-(/ - Z)e/ay - f J)/r 3I } 

wherein 

7 2 (/) : formula defining the previous time history of the intensity 
of the laser pulse 

Delay : time delay of the starting point of the upward edge rela- 
tive to a zero point in time 

t e : difference in time between the beginning of the down- 
ward edge and the beginning of the upward edge 

r 31 : first time constant of the downward edge 

r 32 second time constant of the downward edge 

I 2 (t) is an optional term defining the previous time history of the inten- 
sity of the laser pulse. Preferably, one of the formulae presented above 
should be applied here. 

A further embodiment of the invention provides for the analyzer unit (13) 
to be designed so that it uses approximations of the optical pulse when 
performing a convolution with model functions E{t) for heat transfer for 
the determination of the heat transfer function R(t) : 

R(t)=\l(s)E(t-s)ds 

0 

wherein 

I(t) : pulse profile of the laser pulse 

E(t) : model function for heat transfer 
Then, the thermal conductivity can be determined by deriving it from the 
heat transfer function R(t) . For defining the function of the pulse pro- 
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file /(/) of the laser pulse one or several ones of the aforedescribed 
approximations are applied. 
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In another expedient embodiment of the invention, the analyzer unit 
(13) is configured for performing a convolution of the optical pulse with 
a model function using the Cowan approximation. The Cowan approxi- 
mation, published in the Journal of Applied Physics 34 (1963), p. 926 - 
with the essence of this publication constituting part of this document - 
is based on the following heat transfer equation E(t) as a function of 
time t: 

f kcc 
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E{t) « exp 

wherein 

k 
a 
I 



l + 2j(-irexpf^X^ 



I 2 



J _ 



coefficient descriptive of the thermal losses 
thermal conductivity 
thickness of the sample 



The approximation indicated here starts initially out from a rectangular 
shape of the laser pulse. In combination with the aforedescribed ap- 
proximations for the shape of the laser pulse, an improved precision 
20 can be achieved. 



A further configuration of the analyzer unit (13) consists in its design for 
translucent materials, with application of a heat transfer solution pro- 
posed by Mehling et al. in: International Journal of Thermophysics 19 
25 (1 998) p. 941 - with the essence of that publication constituting part of 
this document: 

E(t) = 2T m ±l Pn exp W0 ^#^^ 
tf " l 2 (j3 2 „+h 2 ) + 2lh 

wherein 
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I 

fin 



sample temperature in thermal balance 

thickness of the sample 

the positive roots of the equation 

radiation losses of the sample with 



a 

£ 



wherein 

Boltzmann's constant 
emission constant of the sample 
ambient temperature 

constant of thermal conductivity 
At temperature differences between the face and the rear side of the 
sample, which are smaller than the sample temperature, the following 
values apply for h: 

{h for even n 

h(l + 2rj) for odd n 

wherein 

7] parameter characterizing the transparency of the sample 

With additional modeling of the thermal losses, in analogy with Cowan's 
proposal, the following equation is achieved: 



20 E{t) = 27; exp(-kat/l 2 y£lp n exp(-aj3 n t)- 



l 2 (tf+h 2 ) + 2lh 
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A further embodiment of the analyzer unit (13) consists in a configura- 
tion of the analyzer unit with application of a heat transfer E(r,t) solution 
published by Cape and Lehmann in: Journal of Applied Physics 34 
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(1963) p. 1909 - with the essence of that publication constituting part of 
this document. In that solution, not only the thermal losses via the face 
and rear side (facial thermal losses) but also the thermal losses via the 
cylinder wall are considered (radial heat losses). 



m=0 



/=0 



^0 = EC m X m £z) I .(r 5 y r )//( S )exp 
wherein 



an 1 



<t>im(t-s) 



V 



ds 



{-\T(2ald)X m 



K+2Y X+ Yt 



D t (r,Y r ) = 



2Y r J 0 (z r (r/r 0 )) 



r 



y;+z! j 0 ( Zi ) 

(l/x) 2 (X 2 Jd 2 +z?/r 0 2 ) 

roots of the equation (X 2 m -Y 2 )tanX m = 2XJ X 

roots of the equation Y r - J Q {z.) = z. • J x (z.) 

facial thermal loss 

radial thermal loss 

Bessel function of the order k 

distance from the centre 

diameter of the sample 



Another embodiment of the analyzer unit (13) consists in a design of the 
analyzer unit for the application of a model function of the heat transfer 
in the case of multiple layers, preferably double or triple layers, respec- 
tively. Here, the thermal conductivity of one layer may be unknown. 

A further embodiment of the analyzer unit (13) consists in a design of 
the analyzer unit for the application of a model function for the heat 
transfer with thermal resistance in the case of multiple layers, preferably 
double or triple layers, respectively. 
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A further embodiment of the invention provides for an analyzer unit that 
is designed to perform a convolution by way of subtraction of the spec- 
trally transformed measured values and subsequent re-transformation. 
Appropriate spectral transformation techniques are the fast Fourier 
transform or also the Laplace transform. 

Description of the Drawings 

In the following, the invention will be described by exemplary embodi- 
ments, without any limitation of the general inventive idea, with refer- 
ence to the drawings. 

Fig. 1 is a general schematic view of an inventive device. 

Fig. 2 illustrates a typical laser pulse of major length. 

Fig. 3 shows a typical laser pulse of medium length. 

Fig. 4 is an illustration of a typical laser pulse of short length. 

Fig. 1 shows an exemplary embodiment of the inventive device. A sam- 
ple (1) is heated to a preselected temperature in a furnace (2). To 
measure the thermal conductivity at this temperature, the sample is ex- 
posed to the energy of a highly efficient light source light (3) in the form 
of short pulses. The energy absorbed by the sample results in heating 
the latter. In correspondence therewith, the infrared radiation (4) emitted 
by the sample changes. This infrared radiation is detected by means of 
an infrared sensor (6) over its time history. The signal of the infrared 
sensor (6) is supplied to an analyzer unit (13) by means of a signal am- 
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plifier (12). This analyzer unit is so designed that it derives from the time 
history of the sensor signal the thermal conductivity of the sample. For 
adjusting the measuring result by elimination of falsifications induced by 
the time history of the optical pulse, a device is provided for measuring 
5 the time history of the optical pulse. This device comprises a decoupler 
element (7) such as a semi-transparent mirror that decouples a prede- 
termined fraction of the radiation from the light source (5) as reference 
radiation (8) and supplies the reference to sensor (9) for measuring the 
laser radiation. This measured value is then supplied to the analyzer 
10 unit (1 3) for determination of the time history of the optical pulse. The 
analyzer unit corrects the signal detected by the infrared sensor (6) by 
performing a convolution with the signal detected by the sensor (9) for 
measuring the laser radiation. A display unit (14) is provided for display- 
ing the result. 

15 

Fig. 2 shows a typical laser pulse of major length. The time unit ms is 
plotted along the horizontal axis. The vertical axis reflects the intensity 
of the pulse in standardized units. The illustrated laser pulse has a rated 
pulse length of 1.0 ms. In such a pulse, three sections can be distin- 
20 guished. It presents 

(a) a high-speed rising starting section 

(b) a generally slightly downward sloping plateau, and 

(c) a high-speed downward slope after cut-off of the pumping light. 
These sections are now approximated by individual functions. 

25 

Fig. 3 illustrates a typical laser pulse of medium length. The time is plot- 
ted in ms along the horizontal axis. The vertical axis reflects the in- 
tensity of the pulse in standardized units. The illustrated laser pulse 
presents a rated pulse length of 0.5 ms. In such a pulse, only two sec- 
3 0 tions (a) and (c) can still be recognized. This pulse can by expediently 
approximated in two segments. 
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Fig. 4 shows a typical laser pulse of short length. The time is plotted in 
ms units along the horizontal axis. The vertical axis reflects the intensity 
of the pulse in standardized units. The illustrated laser pulse has a rated 
5 pulse length of 0.2 ms. In such a pulse, one can only distinguish still two 
sections (a) and (c). This pulse can be expediently approximated in two 
segments. 
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List of reference numerals 



1 


sample 


2 


furnace 


3 


radiation for heating the sample 


4 


infrared radiation emitted by the sample 


5 


source of radiation 


6 


infrared sensor 


7 


decoupler element 


8 


reference radiation 


9 


sensor for measuring the reference radiation 


10 


Laser power unit 


11 


furnace power unit 


12 


signal amplifier 


13 


analyzer unit 


14 


display unit 



